Feedback in the cores of clusters Abell 3581, 
2A 0335+096, and Sersic 159-03 



Of 



Catherine L. Farage 1 , Peter J. McGregor 1 , Michael A. Dopita 123 
cf arage@mso . aim . edu . au ; pet er@mso . aim . edu . au ; madOmso . anu . edu . au 

ABSTRACT 

The cores of massive galaxy clusters, where hot gas is cooling rapidly, appear to 
undergo cycles of self-regulating energy feedback, in which AGN outbursts in the central 
galaxies episodically provide sufficient heating to offset much of the gas cooling. We 
use deep integral-field spectroscopy to study the optical line emission from the extended 
nebulae of three nearby brightest cluster galaxies and investigate how they are related to 
the processes of heating and cooling in the cluster cores. Two of these systems, A 3581 
and Sersic 159-03, appear to be experiencing phases of feedback that are dominated 
by the activity and output of a central AGN. A 3581 shows evidence for significant 
interaction between the radio outflows and the optical nebula, in addition to accretion 
flows into the nucleus of the galaxy. X-ray and radio data show that Sersic 159-03 is 
dominated by the feedback of energy from the central AGN, but the kinematics of the 
optical nebula are consistent with infall or outflow of material along its bright filaments. 
The third system, 2A 0335+096, is dominated by mass accretion and cooling, and so 
we suggest that it is in an accumulation phase of the feedback cycle. The outer nebula 
forms a disk- like structure, ~ 14kpc in radius, that rotates about the central galaxy 
with a velocity amplitude of ~ 200 kms -1 . Overall, our data are consistent with ongoing 
AGN-driven feedback cycles occurring in these systems. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: clusters: individual 
(A 3581, 2A 0335+096, Sersic 159-03) — galaxies: ISM — techniques: imaging spec- 
troscopy 



1. Introduction 

Extended emission-line regions are a character i stic feature of the brightest cluster galaxie s 



(BCGs) in many galaxy clusters (e.g.. lHeckmanlll98ll : ICowie et al.lll983l : ICrawford fe Fabianlll992l ). 
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These nebulae often have filamentary structures and may reach to tens of kiloparsecs from the 
galaxy center. Their spectra feature strong low-ionization emission lines, similar to those from the 
Low- Ionization Nuclear Emission - line Regions (LINERs) that are observed in many galaxy nuclei 
(|Heckmanlll98ol : iFirippenkolhood ). 



BCGs are giant elliptical galaxies, and commonly cD galaxies with extended stellar envelopes. 
They are the most luminous and most massive of the cluster galaxies and are located at, or close to, 
the gravitational center of the host cluster. This position contributes to a rich history of accretion 
and mergers, and also means they are permeated by the densest, X-ray-emitting gas of the ICM 
and may be the destination of any material cooling from this medium. 

In cooling-flow, or cool-core, clusters the X-ray surface brightness profile has a strong central 
peak, indicating a core of high density and reduced temperature. The cooling time in this core is 
short relative to the lifetime of the cluster, so there is the potential for massiv e cooling flows to 
have formed ([Fabian fc Nulsenlll977l ; ICowie Binnevlll977l ; iPeterson et al.ll2003l ). However, X-ray 
spectroscopy has shown that most of the gas cooling is being offset by some source of heating in the 
core. This source is not certainly associated with the activity of the central galaxy, but episodic 
outbursts of radiative and mechanical energy fro m active galactic nuclei (AGN) in the BCGs are 
currently the favored heating mechanism (e.g., see iBegelmanl 12004 : iMcNamara fe Nulsenl 120071 ). 



The line-emitting filaments that surround many BCGs may be symptoms of the interactions 
betw een the cooling cluster gas a nd the sources of mass and energy feedback in the central galaxy 



(e.g., IMcNamara Nulsenl 120071 ) . As such, cluster cores are effective places to study feedback in 



nearby galaxies, through the excitation, kinematics, and origins of the nebulae. 

Observations reveal clear links between the extended emission-line regions and their environ- 
ment. BCGs are more likely to host extended emis sion-line nebulae when they are situated clos e to 
the X-ray center of a strongly cooling cluster core ([Edwards et alJl2007l ; ICavagnolo et al.ll2008l ). In 
addition, BCGs that host more luminous emission-lin e systems tend to be found in more lum inous 



X-ray cluster cores with higher central cooling rates (Edge 



2001 



McDonald et all l20ld . 1201 lh , are 



generally a ssociated with larger masses of cold gas ([Edge et al.ll2002l ). and have higher rates of star 
formation ([Crawford et alJIlQQgl i lO'Dea et al.ll2008l ). This suggests that the presence, or excitation, 
of the extended nebulae in cooling-flow BCGs is directly related to the cooling gas in the cluster 
core. 

Observational results also highlight a connection between AGN output, star formation, and gas 
cooling in BCG environments. For example, there appears to be a minimum threshold in central 



X-ray cooling rate that is re c 
BCGs l|Cavagnolo et al.ll2008l : 



uired for the onset of both nuclear activity and star formation in 



Raffertv et al J 12008 ) 



Models of episodic AGN heating and feedback provide at least a qualitative explanation for how 
the h e ating and regulati o n of cooling flows may occur (e.g., Binnev fc Taborll995 ; Ciotti Sz Ostriker 



20011 : I.Tones et al.1 1 19971 : iTucker fc Davidl 119971 : iPizzolato & Sokerl [2005). In this paradigm the 



BCG and its environment experiences cyclic episodes of cooling and inflow that feed nuclear 
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act ivity, and heating and o utflow powered by the AGN. In the cold feedback model described 



by iPizzolato Sokerl (|2005l ) , feedback occurs throughout the cool inner region of a cluster core 
(over radii of ~ 5 — 30kpc). An outburst from the AGN in the central galaxy ejects a large 
amount of energy into the surrounding medium. Bipolar jets form and drive expanding plumes 
of hot, radio-emitting plasma that rise buoyantly in the cluster atmosphere. These interact with 
the surrounding galaxy and environment, and displace and heat the cooling cluster gas (e.g., see 



McNamara &: Nulsenl 120071 . for a review of AGN heating mechanisms). This heating slows fuel 
supply to the AGN and as outflow activity decreases, material in the cluster core begins to cool. 
Clumps that become sufficiently dense fall inwards and, with some initial angular momentum, may 
form an accretion disk. This produces a reservoir of cold gas in the galaxy; some of the cooling 
material contributes to star formation, and some is accreted into the black hole. A large mass of 
accreting material will trigger another burst of nuclear activity and continue the feedback cycle. 



Wilman et al.l (|2009l ) suggest that the diversity of properties of the multi-wavelength emission 
from cluster cores might be explained, at least in part, because they are observed at different 
phases of such feedback cycles. They discuss several BCG systems and propose that one is in an 
early cooling phase: having gas motions in the filaments consistent with free-fall into the galaxy, 
weak radio emission and a strong X-ray cooling rate. They place two others, hosting stronger and 
more extensive radio sources, in post-AGN outburst phases. Such studies of the multi-phase gas 
surrounding BCGs, particularly including integral-field spectroscopic data that provide kinematic 
information about the emission-line regions, will help us to understand which, if any, properties 
of the extended emission systems are evidence of feedback activity. Further, this may allow us 
to definitively place individual BCG systems within the sequence of a heating/cooling cycle, and 
potentially to constrain the physics of feedback models in clusters. 

We have previously presented a deep, integral-field spectroscopic study of NGC 4696, the 
central galaxy of the Centauru s cluster, that prov ided a detailed picture of the kinematics of the 



central emission-line filaments (jFarage et al.ll2010i ). This system may be in a period of accretion 



following the infall of a small cluster galaxy. We found that models of ~200 kms -1 shocks, 
propagating through pre-ionized, dense clouds, can reproduce the optical emission-line spectrum. 
The results also emphasize that mergers and galaxy interactions are an important component of 
the activity in cooling-flow BCG nebulae. 

Here, we study a small sample of other low- and intermediate-luminosity cooling-flow BCG 
filament systems. These are the first optical integral-field spectroscopy of the nebulae surrounding 
the BCG of A 3581 and Sersic 159-03, and the first to cover the full extent of the emission-line 
regions in that of 2A 0335+096. We investigate the kinematic and excitation properties of the 
optical emission surrounding the galaxies and combine our results with those of studies at X-ray 
and radio frequencies to make inferences about the state of feedback activity in these systems. 

Throughout this study we adopt the cosmological parameters Hq = 71 kms -1 Mpc -1 , Qm = 
0.27 and f^A = 0.73, based on the five-year Wilkinson Microwave Anisotropy Probe ( WMAP) results 
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(|Hinshaw et al.1 12003 ) 



2. Brightest cluster galaxy sample 

We have observed three brightest cluster galaxies that are known to host extended nebulae 
that pr oduce strong low- ionization optical emission lines; e.g., [N II], [O II], [S II], [O I], and 
[O III] (j Johnstone et al.lll987l ; iRomanishin & Hintzen] Il988l ; brawford et al]|l999l ). They are all 
well-known examples of such systems. 

The galaxies represent the low-to-intermediate end of the range of total emission-line luminosity 
of BCGs. The ROSAT Brightest Cluste r Sample (an X-ray-selected sample of BCGs from the 
ROSAT All-Sk y Survey; Allen et al. 19921 ) covers a range of integrated Ha luminosities from 10 39 
to 10 43 ergs" 1 ([Crawford et al.lll999h . The BCGs in A 3581, 2A 0335+096, and Sersic 159-03 have 
Ha luminosities of ~ (5 — 50) x 10 40 ergs -1 . Table [J compares some observed properties of these 



galaxies, together with NGC 4696 (iFarage et all 12010 ). Throughout the rest of this paper we refer 



to the central galaxies by the host cluster names listed in this table. 



2.1. Abell 3581 



A 3581 is a relatively nearby (z ~ 0.022) cluster, however the extended nebula that sur r ound s 
the BCG has been little studie d since early imaging and spectroscopy by iRobinson et al.l (jl987l ) 
and iDanziger &: Focardil (j 19881 ) . The radio and X-ray emission in the cluster core have been in- 
vesti gated in some detail, though, using observations from the VLA, XMM-Newton, and Ch andra 



(e.g., Ijohnstone et al.lll998l : iDunn k. Fabianl 12004 : 1 Johnstone et alJl2005l ; ISanders et al.ll2010l ) 



The BCG is a strong radio source, identified as PKS 1404-267. Two irregular lobes of 1.4 GHz 
emission extend ~10" (4.6 kpc) to the east and west of a compact central source that coincides with 
the optical centroid of the galaxy (j Johnstone et al.lll998l ). Both radio lobes appear to have suffered 
some disruption: there are distortions or bends in the orientation of the extended structures. 

The 1.4 GHz radio lobe s coincide with two cavities in the hot X-ray halo gas of the cluster 
core (j Johnstone et al.ll2005l ). Such structures are seen in many BCG cooling- flow systems and 



are assumed to be evidence that the radio lobe s displace the surrounding h ot gas as they expand 
and propagate away from the central galaxy (see iMcNamara fc Nulsenll2007l ) . iRaffertv et al.l ((2006) 
calculated the total power required to inflate the cavities seen in the X-ray-emitting ICM of A 3581. 
They found that the cavity power and cooling luminosities approximately balance in this system. 
The ratio of cavity power to the cooling X-ray luminosity in the core is: Pcav/^ICM ~ O.8I04. 
X-ray spectroscop y indicates a relatively low rate of gas coo ling to low temperatures in the core: 
M < 10M o yr _1 (jVoigt fc Fabianl 12004 : ISanders et al.ll2010l ). The X-ray and radio observations 
suggest that most of the ICM cooling is offset by the current AGN feedback activity in this system. 
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2.2. 2A 0335+096 



2A0335+09qJ is a compact cluster with a luminous X-ray core ( Sanders et al.1 12009 ). The 
BCG has cD morphology, hosts a weak radio source, and is surround e d by an extensive optical 
emission-line nebula that was first imaged by iRomanishin &: Hintzenl (|1988l ). They reported a 
clumpy, 12 kpc-long bar of line emission and a 2 kpc disk-like feature that they characterized as 
possible accretion disk inclined less than 17° from edge on. 



In optical integral-field spectroscopy of the central 7" x 9" of the nebula, lHatch et al.l (|2007l ) 
resolved two knots of emission within the inner bar of Ha emission. Both knots are less than ~ 2" 
(1.4 kpc) from the galaxy core and show distinct line-of-sight velocities separated by ~ 200 kms -1 . 
These stru ctures were als o seen in detailed near-infrared integral-field observations of the central 
3" x 3" by IWilman et al.l (|201ll ) . A more diffuse region of emission extends from the northwest 
knot towards a nearby galaxy, ~ 6" (4 kpc) away on the sky to the northeast of the nucleus. The 
line-of-sight velocity increases smoothly from ~ —250 kms -1 app roximately 1" sout h of the BCG 
nucleus to ~ 200 kms" 1 at the position of the companion galaxy. Hatch et al. (j2007l ) suggest that 
the secondary galaxy may have disturbed a reservoir of molecular gas in the BCG to produce this 
extension of the inner nebula towards the northwest. 



Long-slit spectroscopy by iGeldermanl (j 19961 ) and iDonahue et al.l (|2007l ) provided evidence that 
the large-scale bar of line emission that extends over more than 6 kpc to the northwest and southeast 
is kinematically distinct from the inner emission structure. The line-of-sight velocity difference 
between the outer regions of the large-scale bar is reversed with r espect to the velocity gradient 
across the inner disk between the galaxies. IDonahue et al.l (|2007l ) noted that a two-dimensional 
velocity map is needed to characterise the velocity field of the nebula. 

In 1.4 GHz and 5 GHz radio images of the system, two diffuse, amorphous lobes of radio 
emission extend ~ 10 — 15" (~ 7 — 10 kpc) from galaxy ce n ter along a p o sition angle of 60° 



(jSarazin et al.l Il995l : IDonahue et al.l 120071 : ISanders et al.1 |2009j) . ISarazin et al.l (|1995l ) suggest that 
these lobes are associa ted with oppos i tely-d irected radio jets and have been disrupted by the high- 
pressure ambient gas. Sarazin et al. (j 19951 ) measure a spectral index of a ~ —0.3 on the galaxy 
nucleus, and a ~ — 1 in the two radio lobes. Based on steepening as a result of synchro tron losses 
(from a ~ —0.5), and assuming an equipartition magnetic field, IDonahue et al.l (|2007l ) estimated 
an age of 25 — 50Myr for the radio lobes. They proposed that an interaction with the nearby 
companion elliptical galaxy may have triggered the current episode of radio activity. 

There is evidence for an active history of interaction between the AGN and ICM in the core of 
2A 0335+096. Studies with XMM-Newton and Chandra have identified a complex system of X-ray- 
emitting structures in the core, including filaments, cool clumps, a metal-rich spi ral, and at least 



five distinct cavities at varying distances and position angles relative to the BCG (jMazzotta et al 



1 identified as RX J0338. 7+0958 in some studies 
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20031 : IBirzan et alJ 12004 : iKaastra et al J 12004 : ISanders & Fabianl 12009 : ISanders et al J 120091 ). The 



radio lobes are associated with weak cavities in the X-ray gas, but the most prominent cavities 
coincide with very extended, diffuse lobes of 1.4 GHz emission approximately 20 — 30kpc from 
the galaxy core on the sky. The complex sys tem of cavities presu mably originates from multiple 
previous generations of AGN radio outbursts (panders et alj |2009). 

X-ray observations indicate a relativel y high mass cooling rate in the core: 30 — 8OM0yr 



(Raffertv et al. 


2006: 


Voiet & Fabianl 


2006 


), and 


power output 


Raffertv et al. 


2006; 


Birzan et al. 



emitt ing material is correlated closely with the large-scale bar of optical line emission ([ Sanders et al 



2009), so there is cooling X-ray gas associated with extended nebula. This system currently appears 



to be dominated by cooling rather than AGN activity. 



2.3. Sersic 159-03 



bersic 159-050 is a relatively low-mass cool-core cluster that has been the subject of many 
observational studies. The core of the cluster has a complex structure at all wavelengths and shows 



signs of strong AGN feedback (jWerner et al 



2011 



Birzan et alJl2008l ). 



Previous imaging ( I 



2010 



Oonk et al 



2010 



ansen ct 



31 



2000 



McDonald et al 



Jaffe et a 



201(1 



2011 



2005; 



Crawford &: 



Werner et al 



2011 



'abianll99 2: McDonald el al 



in the optical and near- 
infrared shows several bright filaments surrounding the BCG. A curved filament extends north of 
the BCG nucleus to a radius of ~ 35 kpc, and comprises a number of filaments and clouds. Hubble 
Space Telescope (HST) observations of the galaxy revealed two lanes of dust absorption that are 



associated with the emission- line structures (jWerner et al.ll201ll ). Faint ultraviolet emission is also 



extended along the filament. The presence of ultraviolet radiation and dust may be evidence for 
star formation associated with this structure. 

Sersic 159-03 is associated with a complex rad i o source. Two inne r lobes of radio emission 



are seen at 5 GHz and 8.5 GHz (Bir zan et al. 



2008 



Werner et al 



2011 



The 5 GHz emission is 

centrally concentrated and the relatively faint lobes are irregular and asymmetric. However, the 
radio emission is considerably more extended at lo wer frequencies. Emission at 300 MHz is detected 
at radii of over an arcminute (jBirzan et al.ll2008l ). The radio source has a very steep power-law 
spectrum (a ~ —2.2) in the extended emission to th e east, northeast, and northwest of th e core, 



Birzan et all 12004* : IWerner et alJboilh . The 



and is more shallow in the nucleus (a ~ —1 to —1.5; 
most recent radio out put (assumed to be the high-frequency emission) is extended to a radius of 
~ 10 kpc in each lobe. IBirzan et al.1 (J200J) suggest that the radio source has either restarted, or is 
fading away. 



2 also designated Abell S1101 (A S1101) 
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Birzan et al.l (|2008l ) also identified indistinct, low-contrast X-ray cavities that lie beyond the 
emission-line filaments and are associated with diffuse, low-frequency radio emission. The power 
associated with inflating these cavities is high, approximately three to four times the luminosity of 
the cooling X-ray gas in the cluster core (jRaffertv et al.l l2008). X-ray spectros copy shows that the 



ICM in the core is cooling to low temperatures at a rate of only < 30 Mp, yr 1 ( Werner et al.ll2011 



McDonald et al. 



2010 



Rafferty et al. 



200ft IVoigt fc Fabianl 12004 ). 



In a recent multi-wavelength study of the BCG core, IWerner et al.l ([2010] ) found evidence for 
strong interaction between the AGN and hot gas of the cluster. They suggested that AGN feedback 
activity has removed most of the cooling X-ray gas from the core (within r < 7.5 kpc). They also 
proposed that the filamentary nebula has been uplifted from the central galaxy by the radio lobes, 
but that gas in the filaments is now largely removed from the impact of the radio jet and is cooling 
and forming stars. 



Observations and data reduction 



The galaxies were observed with the Wide-Field Spectrograph (WiFeS; iDopita et al.1 12007. 



2010) and the Australian National University (ANU) 2.3 m telescope at Siding Spring Observatory. 
WiFeS is a double-beam, image-slicing, integral-field spectrograph that records optical spectra over 
a contiguous 25" x 38" field of view. This spatial field is divided into twenty-five l"-wide slitlets, 
sampled by 0.5" pixels along the 38" length. 

The instrument is a dual-channel device with separate gratings and cameras for the blue and 
red channels. Each camera is equipped with 4kx4k pixel CCD detectors. The low-resolution grating 
configuration, employing the R3000 and B3000 gratings and the RT560 dichroic, was used for the 
observations described here. This provides a continuous wavelength coverage from 3290 to 9330 A 
in the reduced data. The spectral sampling is 0.81 A per pixel in the blue data cubes and 1.3 A 
per pixel in the red data cubes. The typical velocity resolution achieved is ~ 100 kms" 1 FWHM. 

The pattern of science observations consisted of two pairs of 30 minute exposures on a galaxy 
separated by a 30 minute exposure of the sky. The scaled sky frame was subtracted from each of 
the four galaxy frames in this sequence. A dither pattern of 2" offsets was applied to the position 
of each galaxy frame. A nearby reference star was used to autoguide the galaxy exposures. Table [2] 
shows the total integration times obtained for each galaxy and other details of the observations. 

Each set of four galaxy frames and one sky frame was preceded and followed by an observation 
of a spectrophotometric standard star and these were used to calibrate the absolute flux scale of the 
observations. Dwarf stars with spectra that resemble a featureless blackbody were also measured 
during the night and used to remove atmospheric absorption features from the flux standard star 
and galaxy spectra. Cu-Ar and Ne-Ar arc calibration lamp frames, for wavelength calibration, were 
observed regularly during each night, and quartz-iodine lamp flat-field frames were observed at the 
beginning of each night. Twilight sky flat-field frames were obtained on most nights, so that those 
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used in processing an image were never obtained more than two nights from that observation. 



The data were reduced using IRAF data-reduction scripts for WiFeS (jDopita et al.ll2010l ) that 



are based on and directly employ scripts from the Gemini IRAF packag e^. The process of data 



reduction applied to these data is similar to that described in lFarage et al.1 ( 201Q ). to which we refer 



the reader for a more detailed description, but differs in several points that are outlined below. 

A slightly different method was employed to remove the detector bias level from the images. 
Subtracting the bias signal is complicated both by detectors that are read out through four separate 
amplifiers, producing a different bias level in the four quadrants of each image, and by a residual 
shape that is present in the bias level across the detectors. These were removed using polynomial 
fits to the bias signal across unexposed rows on each amplifier in addition to the overscan columns. 

A second-order correction for extinction in the Earth's atmosphere is also included in the 
flux calibration for these data. We estimate that the resulting absolute flux scale is accurate to 
approximately 11%. 

Also, for these data a correction for differential atmospheric refraction was made in software 
after constructing the data cubes, using routines from the Starlink positional astronomy library 
SLALIB in IRAF. 

The average rms error in the wavelength calibration is ~ 0.2 A for the red cubes and ~ 0.1 A 
for the blue cubes. Line-of-sight velocities are presented relative to the local standard of rest, as 
the corrections for the relative velocities of the sources during the observations were smaller than 
the spectral resolution of the data. 

After forming the individual data cubes for each exposure, the pixels were binned by two in 
the spatial direction along the slitlets, to form square 1" x 1" pixels. This improves the signal-to- 
noise ratio and provides a spatial pixel scale that is better suited to the typical seeing during the 
observations (~ 2" - 2.5"). 

Several exposures of A 3581 and 2A 0335+096 were obtained with the field oriented at a 
different position angle (PA) to the majority of the observations. The binned cubes were rotated 
to the intended position angle using the IRAF rotate task, interpolating the pixels linearly. The 
position angles for the combined data cubes are provided in Table [2j 

The individually-reduced science cubes were median-combined after applying spatial offsets 
to align the dithered observations. These offsets were derived from the continuum centroid of the 
galaxy nucleus, measured in an image close to the center of the spectrum in a region that is free of 
line emission. Cosmic rays and bad pixels were also removed by combining the data cubes in this 
way. 



J http : //www . gemini . edu/ sciops/dat a- and-re suits/processing- software 
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3.1. Continuum subtraction 

As seen in Figure [TJ there is a strong contribution to the BCG spectra from the stellar com- 
ponent of the elliptical galaxies. In the central regions of the galaxies, the spectrum of each is 
dominated by this emission, which must be removed to accurately measure the nebular emission 
lines. 

For A 3581, this was done using a spectrum extracted from within the data cube itself, since 
line emission is not observed in all regions spatially across the galaxy. We extracted an integrated 
continuum template spectrum from a position 8" north and 5" west of the continuum peak, within 
a 4"-diameter aperture. 

We could not extract a template galaxy continuum spectrum from within the Sersic 159-03 
and 2A 0335+096 data cubes that was free of line emission and had sufficient continuum signal 
to be useful in removing the galaxy continuum. In these cases, we formed a galaxy continuum 
template spectrum from previous WiFeS observations of the giant elliptical galaxy NGC 4696. An 
integrated spectrum was extracted from a region of high galaxy signal and the line emission was 
removed using the Ga ussian line-profile fi ts to the emission in the continuum-subtracted NGC 4696 



cubes as described in lFarage et al.l (|2010l ). This spectrum was then interpolated onto a redshifted 



wavelength scale for each of the galaxies to form a continuum template spectrum. 

For each of the emission lines of interest, in each spatial pixel the template spectrum was 
scaled to match the flux level of a continuum region close to the line and then subtracted. Figure [2] 
shows example continuum-subtracted spectra from the red and blue data cubes for each galaxy. 
The spectra are extracted from a 4"-diameter aperture located at the galaxy nucleus, as defined by 
the peak of the continuum emission. 

Employing continuum template spectra from the WiFeS data was the most straightforward 
option for removing the stellar continuum, and provides satisfactory results since the continuum 
shape is similar in each of these galaxies. 



3.2. Emission line profile fitting 

The emission-line features in the continuum-subtracted cubes were fit using the Python NMP 



FIT package, an implementation of Craig Markwardt's MPFIT routines for IDL (|Markwardtll2009l ) . 
which applies the Levenberg-Marquardt method for least-squares minimization of data to a model. 
The line-profile models consist of one or more Gaussian functions plus a linear continuum baseline 
to remove any residual continuum level and slope. We constructed and fit the models separately for 
each of the individual emission lines or set of closely-spaced lines, since the continuum subtraction 
is applied locally for each feature or set of features within the cubes. 



The Ha and [N II] AA6548, 6583 doublet were fit together. We constrained the velocity of 
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Ha and [N II] A6548 to be the same as that of [N II] A6583, and the ratio of the [N II] doublet 
line amplitudes to be 3.0. These lines were also constrained to have equal velocity widths. The 
components of the [O III] AA4959, 5007 doublet were fit together and the 5007/4959 A amplitude 
ratio fixed at a value of 2.88. Both lines were also constrained to have the same relative velocity 
centroid and line width. The model of the [S II] AA6716, 6731 doublet also consisted of two Gaussian 
components that have the same velocities and line widths. The components of the [O II] AA3726, 
3729 doublet are sufficiently broad that they are unresolved in the data cubes and a single Gaussian 
component was fit to this line profile everywhere. The line width is broadened by including the 
two components of the doublet in the single fit. A single Gaussian component was also fit to the 
unresolved [N I] AA5198, 5200 doublet. 

The reported line widths have been corrected for instrumental broadening of 100 kms" 1 
FWHM, which is the mean line width measured from arc lamp frames observed with these data. 
Line widths are given as velocity dispersions, rather than FWHM (full width at half maximum). 
Line-of-sight velocities are measured with respect to the redshift of the BCGs as given in Table HJ 

Extinction corrections were derived from the Ha/H/3 flux ratio for the spectra pixels where Ha 
and H/3 are detected above a surface brightness limit of 2 x 10 -17 ergs -1 cm~ 2 arcsec -2 (equivalent 
to a signal-to-noise ratio of 2 to 3 in the data cubes). We assume an intr insic flux ratio Ha/H/3 
3.1. This value is appropriate for active galaxies, including LINERs (jVeilleux &: Osterbrock 



19871 ). in which there is a collisional-excitation component in the Ha flux. It is slightly higher than 



the classical Case B recombination v alue of 2.86. Integ rated emission-line fluxes are corrected for 
extinction using the reddening law of lOsterbrockl (119891 ) and we adopt a value of R = Ay /E(B — 
V) = 3.1. Where the H/3 emission is faint or undetected, we adopt a constant extinction to correct 
the line fluxes. This is estimated by measuring the integrated Ha and H/3 emission-line fluxes 
over a region of the inner, bright filaments in each galaxy (excluding the nuclear region) where the 
signal-to-noise ratio of both hydrogen lines is greater than five, and obtaining a mean extinction 
value from the integrated flux ratio. The average Ay values estimated and used in this way are 
listed in the third col umn of Table EC Th e estimated extinction contributions from the Galaxy, 
from the dust maps of lSchlegel et al.1 (|1998l ). are also listed in Table [3] for each of the galaxies. 

In A 3581, although the emission-line surface brightnesses are strongest in the nuclear region, 
the line profile shapes complicate the extinction estimate from the Ha and H/3 fluxes. There is 
insufficient signal in the H/3 emission to fit the multiple kinematic components that are seen in 
the strong Ha and [N II] emission lines. In particular, in the galaxy core, a low-amplitude but 
very broad emission-line component provides a good fit to the wings of the strong line profiles, and 
contains a significant amount of flux. We were unable to fit this component in the fainter H/3 line 
profiles, and as a result would overestimate the values of the hydrogen flux ratio in the core. For 
this reason, we apply the mean extinction level from the filaments {Ay = 0.4), as explained above, 
to the fluxes in pixels within a radius of 3" in the A 3581 data cubes. 



Uncertainties in the relative line-flux measurements were estimated using the noise properties 
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of the spectra after continuum subtraction. For a given emission line, the la error in the flux density 
per pixel was taken to be the standard deviation in a line-free region adjacent to the feature. These 
were combined in quadrature to form the uncertainties in the integrated line flux. 



4. Results 



4.1. Integrated spectral properties 

Figure Q] presents spectra integrated over a circular 5"-diameter aperture at the peak of the 
galaxy continuum emission for the three galaxies. In each case the galaxy continuum dominates 
the spectrum and is typical of evolved stellar populations in the giant elliptical galaxies. The three 
spectra show strong emission features that are characteristic of LINER emission. Similar features 
are detected in each spectrum. After removing the stellar continuum component, the following 
emission features are detected in each galaxy: forbidden emission lines [O II] AA3726, 3729 (un- 
resolved), [O III]AA4959, 5007, [N I] AA5198, 5200 (unresolved), [O I] AA6300, 6363, [N II] AA6548, 
6583, and [S II] AA6716, 6731, and Balmer hydrogen recombination lines H7, H/3, and Ha. There 
is also a marginal detection of [Ne III] A3967 in the nuclear spectrum of A 3581, and of He I A5876 
in the outer parts of the nebulae of A 3581 and 2A 0335+096. 



In each galaxy, and as seen in other BCG filament systems (e.g., NGC4696; iFarage et al.ll2010l ). 
the morphologies of all the optical emission lines are similar, though detected over different extents 
depending on the strength of the line. T he kinematics are a lso consistent with being the same in all 
the detected lines (e.g., as in Figure 4 of lFarage et al.ll2010l ). Everywhere that the density-sensitive 
[S II] AA6716, 6731 doublet can be measured, the line flux ratios are consistent with densities of the 
[S Il]-emitting gas in the low-density limit: that is, n e < 200cm -3 . 

The total Ha and [N II] A6583 emission-line fluxes measured from the data cubes and the 
inferred luminosities are listed in Table [3) We can estimate an upper limit of the mass of emitting 
i onize d hydrogen in the filaments from the Ha emission-line luminosity, L(Ha), following lOster brock 
(|l989h : 

L(Ha) m p 



M H+ 



(1) 



n e a|f Q h i^E a 

The derived masses are listed in Table assuming an electron density of n R < 200 cm" 3 . We 
adopt an effective recombination coefficient of a|^, = 1.17 x 10 -13 cm 3 s -1 ( Osterbrock!ll989l ). This 
is appropriate for Case B emission at a temperature of T = 10 4 K. The ionized hydrogen masses 
are lo w compared to the mass of molecular material inferred in the filaments. For example, Edge 
( 200ll ) estimate a molecular gas mass of ~ 4 x 10 9 Mq in 2A 0335+096, and find similarly large 
masses in a number of other cooling-flow BCGs. This is consistent with the line emission arising 
in thin layers of ioniz ed gas surrounding dense, cool molecular c ores that contain most of the mass 



in the filaments (e.g., I Salome et al 



2011 



Johnstone et al.1 120071 ). 
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The integrated [N II] /Ha flux ratios are shown in Table [3] with the total luminosities. We note 
that our data for A 3581, 2 A 0335+096, Ser s ic 159 -03, and NGC 4696 show good agreement with 



the correlation observed by ICrawford et al.l (|1999l ) of decreasing [N II] /Ha flux ratio with total 



Ha luminosity in BCG nebulae. Table 0] presents line flux measurements for the emission features 
detected in the spectra of the three galaxies, from a 5"-diameter aperture centered on the galaxy 
nucleus and also in an off-nuclear position in the filaments. For the latter, a bright region of the 
filaments was selected that is sufficiently far from the nucleus to avoid any significant contribution 
to the emission from the AGN. The high [O IIl]/H/3 ratio in the nuclear spectrum of A 3581 is likely 
to be the result of radiation from the AGN, but, as expected, all of the line ratios are consistent 
with a LINER-like emission spectrum. The optical line ratios in the nebulae and the integrated 
emission properties of these objects are typical of BCG systems. 

The morphology and kinematics of the optical line emission detected in the observations of 
each of the three galaxies are discussed individually in the following sections. 



4.2. Abell 3581 

Figure [3] shows integrated continuum and [N II] A6583 emission-line images of the BCG in 
A 3581 from the WiFeS data. The line emission extends over 15 — 20" (7 — 9kpc) towards the 
east and west of the galaxy nucleus. Within the angular resolution of the observations, the peak 
emission-line surface brightness is coincident with the galaxy nucleus. There is unresolved line 
emission in the galaxy core with a peak surface brightness ~ 2 — 3 times higher than in the 
extended emission regions. 

The right panel of Figure [3] compares the 1.4 GHz radio image of this source with the optical 
line emission map. The line emission has a similar orientation and extent on the sky as the extended 
radio emission that is associated with the galaxy. The radio lobes are surrounded by filaments of 
line emission and the brightest regions of radio emission largely coincide with regions where the 
optical line emission is faint or undetected. Based on the morphology of the line emission, there 
appears to be a close association between the radio lobes and line-emitting gas in this galaxy. 

The kinematics of the nebula are complex. There are steep gradients in the line-of-sight velocity 
distribution. There are also several regions where the Ha and [N II] emission-line profiles are clearly 
composed of multiple kinematic components. In particular, within a radius of approximately 2 — 3" 
of the galaxy core, the line profiles show symmetric blue and red wings. These are fit well by a 
second broad, low-amplitude Gaussian component, with a velocity dispersion of 750 — 800 kms -1 
and centered at the same velocity as the strong, narrow line component. This component of the line 
profiles may represent emission from the broad-line region (BLR) of the active nucleus. Within a 
radius of 3" of the nucleus, the profile fits to the narrow and broad velocity components of the lines 
show a uniform gradient in the line-of-sight velocity across the nucleus. The velocity amplitude is 
approximately ±80 kms -1 in this region. 
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The kinematic features and the overall distribution of velocity in the emission-line regions 
are most clearly illustrated using velocity channel maps. Figure U] presents velocity slices across 
the brightest emission line in the spectrum, [N II] A6583. The panels in this figure are formed by 
selecting bins of constant velocity width and summing the wavelength slices in each bin to form an 
image. Each image shows the spatial distribution of gas in the nebula with line-of-sight velocities 
in a given 100 kms -1 range, with respect to the redshifted wavelength of the line. The 1.4 GHz 
radio map is also overlaid on several of the panels for comparison. 

We measure a large range of line-of-sight velocities in the nebula: Av ~ 900 kms -1 . The 
highest velocities, above ±300 kms™ 1 , occur close to the center of the nebula, within ~ 5" (~ 
2.3 kpc) of the nucleus. There is a cloud of high- velocity blueshifted emission northeast of the 
galaxy core where the line-of-sight velocities peak at > 500 kms -1 . This appears as an additional 
blue-shifted peak or wing in the spectral line profiles from the region. At low velocities, between 
approximately ±100 kms -1 , the emission nebula consists of roughly V-shaped arms oriented close 
to the radio axis. This structure is reminiscent of shells or filaments of gas surrounding bi-polar 
cavities. 

The bright spiral filament to the south of the eastern radio lobe produces redshifted emission 
with velocities increasing toward the nucleus and reaching a peak of ~ 450 kms -1 within several 
arcseconds of the core. This filament encircles the southern edge of the radio lobe, and the knot 
at the peak of the radio emission in this lobe lies just interior to the inner edge of the filament. 
Again, the alignment suggests that the shape of the filament is related to that of the radio lobe. 
One possibility is that the emission-line gas may be uplifted by the expansion of the radio bubbles. 
Since we do not have information about the orientation of the radio lobes or emission nebula, we 
can not distinguish inflow from outflow on the basis of the line-of-sight velocities. However, if 
the motion of the emission-line gas is largely driven by the propagation of the radio outflows, the 
velocities are directed outward from the galaxy center. This is consistent with the mostly redshifted 
emission seen in the eastern region of the nebula and blueshifted velocities of the filaments to the 
west. The highest velocity gas in the vicinity of the western radio lobe occurs in a bright clump 
of emission located where the radio lobe changes direction and extends to the south. This may be 
evidence of an interaction of this type occurring there. 

However, in this scenario it is difficult to explain the cloud of high- velocity blueshifted emission 
to the northeast of the galaxy nucleus. The emission lies along the northern edge of the radio 
source. There is emission-line gas with high, but oppositely-directed, velocities located roughly 
symmetrically to the northeast (blueshifted) and southeast (redshifted) of the nucleus; seen by 
comparing the \v\ = 250 — 350 kms -1 velocity slices in Figure [H This is not consistent with an 
outflow driving the kinematics in (both) these structures. 



It is possible that the material in the spiral-shaped filament to the southeast of the galaxy is 
flowing inward on an orbit-li ke trajectory. This is reminiscent of the motion of gas observed in the 
main filament of NGC 4696 (jFarage et al.ll2010i ). In this case the high- velocity blueshifted emission 
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to the northeast of the nucleus may be gas on the approaching side of this spiral orbit, after it 
passes behind the nucleus. The maximum velocity of ~ 450 kms" 1 occurs in the blueshifted clump 
approximately 3.5" (1.6 kpc) from the position of the nucleus. As suming circular orbita l motions 
implies an enclosed mass within this radius of ~ 5.9 x 1O 1O M0. iRaffertv et al.l ( 20061 ) estimate 



that the bulge mass of this galaxy is (5.7 ± 0.5) x 10 M & so this inte rpretation implies t hat 



approximately 10% of the bulge mass is located within a radius of 1.6 kpc. iGovoni et al.1 (|2000i ) fit 
a de Vaucouleur's profile to the i?-band surface brightness of the galaxy and obtained an effective 
radius of R e = 29" ( 13.2 kp c). A de-projection of this profile assuming spherical symmetry, as 
performed by lYounel (|1976l ). indicates that the fraction of the total luminosity found within a 
radius of 3.5" is ~ 5%. This implies that a similar mass fraction lies within this radius, assuming 
a constant mass-to-light ratio throughout the galaxy. This is similar to the result derived from the 
kinematics and so the measured velocities plausibly arise from material on an orbital trajectory 
and could represent infall. 

In the outer regions of the filaments where the single Gaussian line-profile fits are unambiguous, 
the gas velocity dispersion is relatively uniform in the nebula. Several ridges of broader lines are 
consistent with being the result of beam-smearing, where the line-of-sight velocity gradient is high 
and pixels sample a range of velocities, causing broadening of the line profiles. The average velocity 
dispersion in the outer nebula is ~ 90 rfc 15 kms" 1 . 

The ratio of [N II] to Ha flux is also quite uniform in the outer filaments, with values of 
~ 1.0 — 1.4. From the dual-component Gaussian fits, we find a slight increase in the [N IlJ/Ha 
ratio in the core: the narrow line components indicate a flux ratio of approximately 2.0 in the 
central pixels. The [O IIl]/H/3 line flux ratio shows a distinct increase at the nucleus, from values 
of ~ 0.5 — 1.0 in the spiral filament to ~ 2.5 in the nucleus. Both line flux ratio increases are 
attributable to a contribution to the ionizing spectrum from the radiation field of the AGN. 



4.3. 2A 0335+096 

Figure [5] shows a map of the [N II] A6583 emission detected in our observations of 2 A 0335+096 
and an image of the continuum emission from the red data cube. The neighboring cluster galaxy 
to the northwest of the BCG is indicated in the figure. Our integral-field data map the kine- 
matic pr operties across nearly the full extent of the nebula that has been detected in narrowband 



imaging (jRomanishin &: Hintzenlll988l ). The detected line emission extends over almost the entire 



26 kpc x 26 kpc field of view. 

Figure [6] shows the properties of single Gaussian line fits to the [N II] A6583 line. However, there 
are several regions in the nebula where the lines have more complex profiles than a single Gaussian 
kinematic component. This increases the uncertainties in the line flux, dispersions, and velocity 
centroids in those regions, and the increased velocity dispersions surrounding the BCG nucleus 
and to the north of the companion galaxy are the result of this kinematic structure. The velocity 
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channel maps in Figure [7] also illustrate the line-of-sight kinematic structure in the emission-line 
regions. 

A complex system of filaments and clouds extends beyond the centr al BCG and companion 
galaxy . The brightest emission is confined to the bar structure described by lRomanishin fc Hintzen 
(j 19881 ) , which extends across the field of view is oriented at a position angle of approximately 140° . 
The velocity map in Figure [6] c learly demonstrates the kinematic stru cture that was seen in this bar 
in early spectroscopic studies (|Geldermanlll996l ; iDonahue et al,ll2007l ). and reveals the distribution 
over the entire emission region. We confirm that th e outer parts of the nebula appear to be 
rotating about the center of the BCG, as proposed by IDonahue et al.l (|2007l ). The rotation has a 
peak circular velocity of ~ 200 kms -1 at a radius of r ~ 15" (lOkpc). The velocity gradient across 
the inner region of the nebula between the two galaxies is superimposed on this large-scale velocity 
field and counter-rotates with respect to the outer bar with an amplitude of Av ~ 130 kms" 1 . The 
line-of-sight velocities associated with the relative motion of the BCG and companion galaxy share 
this counter-rotational sense of motion with the inner disk-like structure around the nucleus. 

The right panel of Figure [5] compares the morphology of the 1.4 GHz radio source associated 
with the galaxy to the optical line map from our data. There are faint filaments of line emission 
that extend to the n ortheast and southwest of the galaxy nucleus and encircle the radio lobes. 
Donahue et al.l (|2007l ) highlight two circular regions as being delineated by line emission and co- 
inciding approximately with the radio lobes. These can be seen in the central panel of Figure [5] 
where the contours outline roughly circular depressions in the emission-line surface brightness to 
the northeast and west of the galaxy. The line-of-sight velocities in these structures are low, which 
could indicate that they have only low inclination from the plane of the sky and are oriented close 
to the rotation axis of the emission-line gas. 

The ratio of [N II] /Ha is approximately uniform along the extent of the bar of bright line 
emission, including between the BCG and companion galaxy, with a value of 1.0 ± 0.1. The 
measured ratios of flux in the components of the [S II] AA6716, 6731 doublet are also consistent 
with being uniform over the brightest emission regions where the signal-to-noise ratio is high, and 
indicate electron densities of < 200 cm -1 . 

The [O III] A5007 line is detected only from the inner region of the nebula between the nuclei 
of the two galaxies, but is extended along the axis of the bar and shows a similar velocity structure 
as the brighter Ha and [N II] lines, consistent with arising in the same extended filamentary gas. 
The [O IIl]/H/3 ratio increases at the nuclei of both the BCG and companion galaxy, with values 
of ~ 0.7 ± 0.1, but is low (< 0.2) in the outer regions of the bar, where H/3 is detected and [O III] 
is not. 



-16- 



4.4. Sersic 159-03 

The extent of the optical line emission detected surrounding Sersic 159-03 is illustrated in 
Figure [8l together with an image of the galaxy in continuum emission. We observe the inner 20" 
(21 kpc) of the bright, curved northern filament, which extends beyond the field of view of our data 
cube. The surface brightness in the filament increases toward the galaxy nucleus and the peak line 
surface brightness is offset by ~ 1" to the east of the galaxy continuum peak. 

The right panel of Figure [8] shows the 5 GHz radio map of Sersic 159-03 overlaid on the [N II] 
image from our data. There is no clear correspondence between the morphology of the optical 
line-emitting filaments and the radio emission, to the degree that the angular resolution of both 
observations reveals. However, as illustrated in Figure (9j the northern filament is extended along 
approximately the same position angle (160 — 170°) as the emission detected at radio frequencies 
from 0.3 to 8.5 GHz. The southern radio lobe is not associated with any detected optical emission 
and there are several filaments that are not located close to any radio structures. 

A single Gaussian function provides a good fit to the emission-line profiles everywhere in the 
data cubes; we do not resolve more complex kinematic components. Figure [TUl shows maps of the 
fluxes, velocity centroids, and velocity dispersions derived from fits to the [N II] A6583 line. Velocity 
channel maps across this emission line are shown in Figure [TTJ 

The line-of-sight centroid velocities of the detected emission vary over ~ 330 kms" 1 . The 
velocity centroid map and channel maps show a velocity shear across the nucleus of the galaxy 
that resembles rotation of the gas in the nucleus, with an amplitude of approximately ±60 kms -1 , 
over a radius of ~ 2" (~ 2 kpc). This provides an estimate of ~ 10 9 M for the enclosed mass, 
assumin g circular rota t ion. T his is consistent with the black-hole mass estimate of (4 ± 1) X 10 8 



Mq bv lRaffertv et al.l ((2006), based on the central velocity dispersion. The axis about which the 



implied gas rotation occurs has a position angle on the sky of ~ 160°. 

The nebula shows high central line-of-sight velocities and velocity dispersions, and decreasing 
velocities outward along the filaments. The velocity from the emission in the northern filament 
decreases from ~ 90 kms" 1 within 1 — 2" of the nucleus, to approximately zero at the outer end 
of the filament. The northern filament diverges into two ridges beyond 5 — 10" from the galaxy, 
and the line-of-sight velocities in the southern component are redshifted by ~ 50 kms -1 relative 
to the north region. This is evidence of more complex structures that are not well resolved in our 
data. The kinematics do not allow us to distinguish between infall and outflow in these filaments, 
since the inclinations are not known. However, the line-of-sight velocities of the inner part of the 
northern filament are consistent with a smooth transition between the velocities in the filament 
and the velocity shear across the nucleus. 

The line-of-sight velocity dispersion in the emission-line gas decreases towards the outer regions 
of the filaments, from ~ 100 kms -1 in the bright inner region of the northern filament, to ~ 
60—80 km s _1 at the outermost edges of the nebula. This could be the result of increased turbulence 



-17- 



in the inner nebula, or the superposition of filaments with more complex velocity structure that 
broadens the measured line widths. 

The line flux ratio of [N II] /Ha peak s in the galaxy cor e and is also generally lo wer in the outer 
regions of the nebula, as was observed by iJaffe et al.l (|2005l ) and lOonk et al.l (|2010l ) . The flux-ratio 
values are also lower in the western filament ([N Il]/Ha~ 0.7) than in the filament to the north 
([N Il]/Ha~ 0.8). Figure [TOl shows the [N II] /Ha flux-ratio map of Sersic 159-03 from our data. 



5. Discussion 

5.1. Interactions between the radio lobes, ICM, and optical nebulae 

Interactions between the line-emitting nebulae, the cooling ICM, and radio outflows that heat 
the environment provide clues to how the filaments are linked to feedback processes in cluster cores. 
For example, several BCGs have been observed where the filamentary structures largely avoid the 



Sabra et al. 
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O'Dea et al. 
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2011 


)■ 


Sabra et al. 



(J2000J) suggest that this is evidence either that the radio lobes displace and entrain the nebular 
gas, or that the relatively dense, line-emitting gas confines the radio lobes and channels the outflow 
through the regions of lowest density in the inter-galactic medium. This raises the question of 
whether the principal effect is of the radio jet and lobes influencing the morphology and kinematics 
of the emission-line gas or vice versa. These issues are linked to whether the emission-line gas 
is infalling and primarily a symptom of gas cooling and accretion, or outflowing and related to 
reheating facilitated by the AGN and radio jet. Of our sample, this type of interaction is particularly 
notable in A 3581, though it may have occurred to some degree in the other two systems. 

In A 3581, the kinematics and the morphologies of the ionized gas and the radio emission 
suggest that the two components are interacting. The shape of the optical emission regions, forming 
bipolar structures that surround the radio lobes, suggests a component of outflow in the nebula that 
is generated by the expanding lobes. Given the coincidence between several clumps of line-emitting 
gas and bends in the orientation of the radio lobes, we conclude that there is significant interaction 
between the structures. The gas with low line-of-sight velocities that surrounds the outer sides 
and edges of the radio bubbles has oppositely-directed velocities at the east and west ends of the 
nebula, and so is probably gas that is being displaced and pushed outwards from the galaxy by 
the expansion of the lobes. The distortion of the radio lobes implies that the radio structures are 
also perturbed by the presence of cold, dense gas in their path. We suggest it is most likely that 
the radio lobes and cooling optical nebula mutually affect one another as the lobes expand, since 
there is evidence for both possible inflow and outflow motions in the nebula; the radio lobes are 
driving out gas in the outer filaments, while there may be gas infalling along the inner, bright spiral 
filament. 
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The structure of the radio source, the X-ray cavities, and the morphology and kinematics that 
we have observed in the optical nebula reveal the system to be in a well-developed stage of AGN 
outflow and heating. The radio source is not overly luminous, but there is significant interaction 
between the radio outflows and the surrounding environment. There is evidence for infall and 
accretion occurring from the line-emitting nebulae, which may be remnant flows that are feeding 
the current nuclear activity. It appears that we are observing a phase of feedback activity in which 
AGN heating balances and possibly dominates the cooling of gas from the ICM. 

In 2A 0335+096, there is no indication of a strong interaction between the optical and radio- 
emitting structures. The filaments to the east and west have some spatial association with the 
radio bubbles: the lobes are coincident with low surface-brightness regions of the nebula, and there 
are faint filaments that outline the edges of these cavities. It is possible that some of this emitting 
gas has been entrained by the expansion of the lobes to their current positions. However, most of 
the emission-line gas in the system appears to be unrelated to the weak radio lobes in this galaxy. 
In particular, the large-scale bar of ionized gas that dominates the structure is not associated with 
any detected radio emission. 



Sanders et al.l (|2009i ) estimate that the total enthalpy associated with the five X-ray cavities 



that they detect in 2A 0335+096 is approximately 5 x 10 59 erg. The cooling lu minosity from 



within the 120 kpc cooling radius of the cluster is approximately 3 x 10 ergs (Bir zan et al 



2004 



Rafferty et al.ll2006l ). Therefore, they note that if the average cooling rate is uniform over 



time, the estimated enthalpy of the cavities would be sufficient to balance the ICM cooling for 
only ~ 5 x 10 7 yr. The current radio sou rce is estimated to have an age similar to this timescale 
([Donahue et al.l 120071 ; ISarazin et al.lll995l ). and estimates of t he average time between outbursts in 
cooling- flow clusters are closer to ~ 10 8 yr ( Birzan et al.ll2004l ). Given the evidence for a number of 
episodes of AGN activity that have produced the observed cavities, the heating has clearly occurred 
intermittently over a significantly longer timescale than 5 x 10 7 yr: at least the length of several 
AGN outburst cycles. Thus, for a significant fraction of the time the core must have been dominated 
by cooling, and it currently appears to be in such a phase of strong cooling and accretion activity. 

The structure and kinematics of the optical emission also provide evidence for material forming 
a rotating disk-like configuration in the outer nebula, possibly being accreted from the cooling gas 
of the ICM. There is a large mass of cool material associated with the extended filaments. CO 
emission-line measurements imply that there is ~4x 10 9 M G of molecular material within 30 kpc 



of the BCG (lEdge 



2001 



Edge Fraverll2003h . 



It appears that cooling activity dominates in the core of 2A 0335+096 and the energy output 
from the radio source is insufficient to balance the ICM cooling at the current time. The radio 
outflows appear to have relatively little influence on the emission-line nebula, and we suggest that 
we are observing a system in a strongly cooling phase of feedback. 



The data for Sersic 159-03 are more difficult to interpret. We can not discount an interaction 
between the optical filaments and radio outflows, but do not see clear evidence for it in our data. 
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Werner et al.l (|201ll ) argue that there has been a strong interaction between the AGN outflows and 
optical-line-emitting filaments. They propose that the material in the north and west filaments 
has been uplifted from the central galaxy by an expanding radio outflow, and assisted in reaching 
the current position relative to the galaxy by the southward motion of the BCG through the 
cluster atmosp here. The BCG is located approximately 8"(8.6kpc) south of the cluster X-ray 



emission peak ( Jaffe et al.ll2005l ) and has a redshift of ~ —500 kms -1 relative to the cluster center. 
This interpretation is supported by metallicities and temperatures measured in the clumpy X-ray 
emission aligned with the optical filament that are lower than in the surrounding ICM. In addition, 
they find that the core of the cluster is largely devoid of the densest X-ray emitting ICM within 
~ 8 kpc, which has presumably been evacuated by the AGN outflows. 

However, there is not a strong relationship between the distribution of the optical line emission 
and the radio lobe structures that are observed. The western filament, in particular, is quite separate 
from the radio emission and seems unlikely to have been affected by the radio outflow activity. Our 
data indicate that the kinematics of the northern filament are equally consistent with inflowing 
velocities along the structure as with outflow, and we suggest that the material in the filament 
could be part of an accretion flow into the central galaxy. The kinematics of the material in the 
filament are consistent with a flow that smoothly joins the rotation seen in gas close to the core. 
The filament material may have been uplifted by outflows associated with an earlier episode of 
AGN activity, have since had time to cool, and is now falling into the central galaxy, potentially 
fueling renewed AGN activity represented by the high-frequency inner radio lobes. 

It is clear that the central AGN has had a significant effect on the surrounding ICM. So although 
the kinematics of the optical emission-line regions do not assist in a definitive interpretation of the 
behavior of the material in the nebula, it is likely that we are observing this system in an AGN- 
dominated phase of feedback. 

To summarise, in A 3581 the structure of the nebula is closely related to the radio outflows and 
it seems clear that AGN feedback at least offsets most of the core cluster cooling in the observed 
phase of activity. In 2A 0335+096 there is evidence for minor interactions between the emission- 
line gas and radio lobes, but predominantly the nebula is not affected by the AGN output, and 
accretion and cooling activity dominates. Sersic 159-03 shows evidence for AGN activity that has 
had a strong influence on the cluster core, and may have affected the cooling material in the nebula, 
particularly the northern filaments, in the past. The optical spectroscopic data provide information 
about the state of the warm gas in the cluster core that can reinforce inferences drawn from the 
X-ray and radio data, but do not necessarily provide an indication of the state of feedback in the 
absence of information about the hot cluster gas and/or radio activity. 
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5.2. Core gas rotation and radio-axis alignment 

In all three BCGs there is evidence for rotation in the emission-line gas close to the core of 
the galaxies. Such motion could represent accretion flows in the nuclear regions that provide a 
channel for feeding material to the central black hole. Furthermore, the axis of rotation within a 
few kiloparsecs of the nucleus is, in all three cases, approximately aligned with the outflow axes of 
the radio lobes associated with the galaxies. 

In A 3581 there is a gradient in the line centroid velocities derived from fits to the narrow 
emission lines in the galaxy core. This resembles a velocity field of rotation about an axis oriented 
at PA ~ 115° and may be evidence for an accretion disk in the core. From the 1.4 GHz radio map, 
the orientation of the radio lobes in the inner part of the jets is approximately 95° on the sky. If 
the bright spiral filament in the nebula of A 3581 does consist of material that is currently falling 
into the central galaxy, as suggested in Section T4.21 it may represent a remaining component of the 
accretion flows in the galaxy core that contribute to the current nuclear activity. The orientation 
of the filament is not well resolved close to the center of the galaxy, but the structure bends 
towards the north in the inner regions and the data are consistent with a trajectory that is close 
to perpendicular to the radio axis as it approaches the nucleus. 



In 2A 0335+096, IWilman et al.l (|201ll ) have mapped the velocities of the H2 emission that 
forms a dumbbell-shaped structure within the central 2kpc. There is a velocity shear across this 
emission region that could represent rotation of gas in an edge-on torus or disk. The axis of rotation 
has a position angle on the sky of ~ 55°, which again is close to the orientation of the radio lobes 
(PA ~ 60° ) . Material from this structure could feed the galaxy nucleus and the orientation of the 
radio outflow axis may have been established by the equatorial accretion flows from this disk of gas 
in the galaxy core. 

If the bar-shaped structure and rotational motion of the outer nebula represent gas settling 
into a large-scale disk, it is likely to be a manifestation of cooling and accretion activity on large- 
scales over the cooling region of the cluster core, occurring during a cooling-dominated phase of the 
feedback cycle. The opposite sense of the rotation in the outer part of the nebula compared to the 
motion of gas in the nucleus may indicate that this activity is decoupled from the phase of accretion 
occurring in the core. This decoupling represents a break in the accretion of material into the inner 
part of the nebula that likely results in an interruption of AGN fueling and a therefore a lapse in 
the feedback process. This may be what allows the cluster core to be largely cooling-dominated at 
the observed period. 

In Sersic 159-03 there is a continuous gradient in line-of-sight velocity across the nucleus that 
resembles rotation about an axis with PA ~ 160°. This axis is close to the axis of the extended 



radio lobes at PA ~ 170°. lOonk et al.l (|2010j) also note that there is evidence from their near- 
infrared, integral-field spectroscopy that the H II and H2-emitting gas within several kiloparsecs of 
the nucleus rotates about an axis parallel to the radio jet. 
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These relationships between the radio axes and the accretion flows in the cores suggest that the 
inflows and outflows are closely connected, as expected in the AGN feedback model. The timescales 
of the implied accretion processes provide a measure of the feasibility of these links. We can make 
estimates about the rate of gas accretion into the galaxy cores if we assume that the observed 
motions represent material in accretion flows that feed cooling material from the nebula into the 
nucleus, and that the bright filaments in Sersic 159-03 and A 3581 trace infalling gas feeding these 
flows. The approximate accretion time through the filaments may be estimated from the average 
velocity v of material along the filament and the distance dfy along the structure into the core: 
i acc dfix/v. The mass inflow rate in the filament is then 



Mm 

10 7 M CT J V 100 km s 



-,-1 



( dm 
\kpc 



(2) 



where Mm is the total mass of material in the filament. 



This total filament mass is uncertain, since the detailed structure and composition of the 
nebulae are not known. Our optical line flux measurements indicate that the total mass of ionised 
gas in the filaments of each in these BCGs is less than ~ 10 6 Mq (see Table [3]). However, this 
probably represents only a fraction of the total m ass, si nce most of the m aterial is likely to exist 
in cool, dense phases (e.g., Johnstone et al. 2007 ). From [jaffe et al.l (j2005l ). the approximate mass 
of warm H2 in the filaments of Sersic 159-03 is on the order of 2 x 10 5 Mq. CO detections imply 
a mass of ~ 1 9 Mq of c old molecular gas associated with BCGs of similar luminosity (including 



2A 0335+096; [Edge 



2001 



). Assuming that the majority of the material is located in the cool, dense 
molecular cloud cores in the filaments, the mass in the inner nebula may be 10 8 — 1O 9 M0. A 
significant fraction of this could be located in the bright spiral filaments that dominate the nebulae 
in A 3581 and Sersic 159-03. We adopt a mass of Mgi ~ 10 7 Mq in these filaments, and an average 
infall velocity of v = 300 kms -1 . The distances along the filaments are l&\ ~ 7 — 15kpc for A 3581 
and ~ 15 — 30kpc for Sersic 159-03, since the degree of projection is uncertain. Then the minimum 
mass- inflow rates in the filaments are on the order of several 0.1 Mq yr _1 , with infall timescales of 
2 x 10 7 to 10 8 yr. 

To consider how these inflows might relate to accretion in the nucleus, we compare the infall 
rates with estimates of the black hole accretion rates required to fuel the observ ed radio activity. 



The 1.4 GHz radio power of the radio lobes is 1. 4 x 10 24 WHz 1 i n A 3581 (jjohnstone et al 



19981), and 1.5 x 10 24 WHz"' for S ersic 159-03 (IBirzan et al.1 12008 ). Using the radio scaling 
relation derived by lBirzan et al.1 (120081 ) for the radio luminosity at 1.4 GHz, £i.4CHz : 



logP cav = (0.35 ±0.07) log L1.4GHZ + (1.85 ±0.10), 



(3) 



the power required to inflate the lobes P cav in these systems is on the order of P c: 



10 42 ergs , though there is significant scatter in this relation. The accretion rate M ; 



= 80 x 
required 

for this output is obtained from M acc = P cav / ( ec2 )> where e is the efficiency of the conversion of the 
gravitational binding energy of the accreting material into jet energy. This is the result of largely 
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unknown physics of AGN jet formation, however a value of e ~ 0.1 is commonly adopted (see e.g., 



Raffertv et al.l 120061 ) . T he required nu c lear a ccretion rate is then approximately 0.01 Mq yr 1 for 



both of these galaxies. iRaffertv et al.l (|2006l ) also estimate accretion rates from energy needed to 
create the physical volumes of the observed X-ray cavities in these clusters. They obtain black hole 
accretion rates of ~ 0.003 and ~ O.lMQyr -1 for A 3581 and Sersic 159-03, respectively. 

That the filament mass infall rates exceed the average black hole accretion rates, particularly 
given the uncertainty in the parameters, emphasises that these structures could provide sufficient 
gas flows to power the AGN in these galaxies. However, this is dependent on the efficiency of mass 
transfer through the inner accretion flows or disks and onto the black hole. The results suggest 
that the filament in A 3581, at least, could fuel the central AGN at the current activity level for 
~ 10 7 yr. Given the similarity of the rates derived for Sersic 159-03, however, mass flow through 
the filaments could, at some times, limit the fuel for nuclear activity if they are the principal source 
of material feeding accretion, in this system. The timescales for infall are similar to the estimated 
frequency of AGN outbursts, 10 to 10 8 years, which also emphasises that the filaments may be 
conduits for material inflow during cooling episodes that trigger episodes of activity. In addition, 
it implies that we might expect to simultaneously observe symptoms of inflow and accretion with 
those of outflow and heating in some systems, as we do in these galaxies. 

In summary, there is evidence in each BCG for an equatorial component of gas rotation in 
the galaxy nuclei, where the bipolar radio lobes, or outflows, are oriented approximately along 
the corresponding polar direction. This is consistent with gas accretion in the cores of all three 
galaxies that may be fuelling nuclear activity and the AGN outflows. It corroborates the model of 
AGN-powered heating and feedback in the cluster cores, though observations with higher spatial 
resolution will be required to better constrain the geometry and properties of the flows. 



6. Conclusions 

We have presented optical integral-field spectroscopic observations of the line-emission regions 
surrounding the BCGs in A 3581, 2A 0335+096, and Sersic 159-0 3. Galaxies such as th ese are 



theorized to be participating in a cycle of negative feedback (e.g., iPizzolato Sz Sokerl 120051 ) . This 
takes the form of energy and mass exchange between the central black hole, host galaxy, and intra- 
cluster medium, in which the heating that results from the nuclear activity in the galaxy can offset 
the energy lost to cooling of the thermal gas in the cluster core. The cycle is self-regulating because 
the AGN output is governed by the amount of mass that cools and is eventually accreted by the 
black hole. 

We have investigated the kinematics and distribution of optical line emission in the extended 
filamentary nebulae, and compared these properties with the extent and development of the radio 
emission in the BCG, and with the properties of the ICM. These data provide evidence of the 
stage of feedback activity that is occurring in these systems, in the framework of the cold feedback 
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scenario. 

The data suggest that A 3581 may be in an advanced heating phase of the feedback cycle. The 
radio lobes are well-formed and appear to have some impact on the surrounding gas, and also to 
have themselves been influenced by interactions with the surrounding medium. There is evidence 
that the expanding radio lobes are evacuating material from the vicinity of the galaxy. Cooling 
and inflow still appear to be occurring in the nebula, though, with gas possibly on infalling orbital 
trajectories along the prominent spiral filament in the central nebula. The kinematics of the gas in 
the galaxy core are consistent with rotation in accretion structures in the nucleus. 

2A 0335+096 is likely to be at the opposite end of the feedback cycle, in a phase of strong 
cooling and accretion, following an episode of AGN activity. A large mass of gas is cooling within 
a wide region surrounding the galaxy. In this system most of the cooling material appears to be 
forming a large-scale disk of gas and dust around the central galaxy. There is evidence for rotation 
and accretion flows on kiloparsec scales in the galaxy core. 

The observations of the optical line emission in Sersic 159-03 show complex structures and 
kinematics, but do not provide a conclusive picture of the influences on the filaments. As observed 
in radio and X-ray data, the cluster core has been strongly affected by energy feedback from the 
central AGN in the galaxy, and likely is experiencing significant heating and only moderate cooling 
and accretion. There is evidence for gas rotation on scales of several kiloparsecs that may represent 
accretion flows into the nucleus. 

The data for these galaxies appear to be at least qualitatively consistent with a model of cyclic 
feedback in the cluster cores, so that, as discussed, we can begin to place these systems within the 
sequence of the feedback cycles that are thought to occur in the cores of galaxy clusters. 
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Table 1. Properties of the BCG sample 



Cluster ID 


BCG ID 


z 


D (Mpc) 




-^rad 


Lk c 


Afbulge d 


A 3581 


IC 4374 / PKS 1404-267 


0.0218 


93 


4.7 


5.4 


0.27 


5.7 


2A 0335+096 


2MASX J03384056+0958119 


0.0346 


147 


49 


1.2 


3.4 


18 


Sersic 159-03 


ESO 291-9 


0.0564 


240 


23 


21 


2.2 


11 


A 3526 


NGC 4696 


0.0099 


44 


1.8 


8 


0.28 


11.6 



a Integrated total luminosity of Ha in 10 40 erg s 1 , from this work and iFarage et al ] l|2O10h . assuming 
the distances in column 4. 

b Radio luminosity in 10 40 ergs -1 , integrated between 10MHz and 10GHz, from lBirzan et aL | ll2008l l, 
except the value for A 3581, which is integrated from a power-law fit to the radio flux data available in the 
NASA Extragalactic Database (NED). 

c X-ray lu minosity of the cluste r inside the cooling radius in 10 44 ergs -1 , integrated between 0.001- 
100 keV, from lRaffertv et alj <2006ft . 

Bulge masses in 10 11 Mq , from lRaffert7 et al. I l l2006h , based on the ij-band magnitudes. 



Table 2. Summary of WiFeS observations 



Host cluster 


Scale 


Dates observed (2009) 


Total exp. time 


PA 




(kpc arcsec - ) 




(hr) 


(deg) 


Abell 3581 


0.455 


May 1-3; 23-24 


11.5 


180 


2A 0335+096 


0.675 


Oct 15-21; Nov 13-15 


11.0 


330 


Scrsic 159-03 


1.074 


Jul 20, 23-24; Oct 16, 21 


11.5 


188 



Table 3. Total Ha and [N II] fluxes and luminosities 



Cluster ID 


^V,gal a 


4„ b 

^V^avg 


F(Ha) 


L(Ha) 


M(H II) 


F([N II]) 


L([N II]) 


[N II] /Ha 


A 3581 


0.20 


0.4 


4.5 + 0.5 


4.7 + 0.5 


< 5.5 x 10 5 


6.8 + 0.7 


7.0 + 0.8 


1.5 ± 0.2 


2A 0335+096 


1.4 


1.5 


19 + 2 


49 + 6 


< 5.8 x 10 6 


18 + 2 


47 + 6 


1.0 + 0.2 


Sersic 159-03 


0.04 


0.5 


3.3 + 0.4 


23 + 3 


< 5.7 x 10 6 


2.9 + 0.3 


20 + 2 


0.9 ± 0.1 



Line fluxes are in units of 10 14 ergs 1 cm 2 , and luminosities in 10 40 ergs 1 . Column 5 shows the estimated 
upper limit on the mass of ionized hydrogen in Mq, for an electron density n e < 200 cm 3 . 

a Foreground Galactic extinction (in magnitudes) from lSchlegel et al. I GUI). 

b Extinction (in magnitudes) assumed in pixels where it can not be estimated from the Ha and H/3 fluxes and the 
assumed extinction law. 
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Table 4. Extinction-corrected line fluxes in selected 5"-diameter apertures 





A 3581 


2A0335±096 


Sersic 159-03 




Nucleus (0,0) 


Filament (-4,-4) a 


Nucleus (0,0) 


Filament (l,17) a 


Nucleus (0,0) 


Filament (-l,4) a 


F(H/3) b 


3.8 ±0.4 


1.5 ±0.2 


8± 1 


4.1 ±0.5 


2.3 ±0.3 


2.2 ±0.3 


Emission-line strength relative to _F(H/3) 


[O II] A3726+3729 


4.9 ±0.2 


6.6 ±0.2 


6.1 ±0.2 


6.0 ±0.4 


6.3 ±0.2 


5.0 ±0.1 


[O III] A5007 


2.16 ±0.07 


0.77 ±0.03 


< 0.63 


< 0.58 


< 0.54 


< 0.37 


[N I] A5 198+5200 


< 0.33 


< 0.40 


< 0.29 


< 0.54 


< 0.31 


< 0.37 


[O I] A6300 


0.66 ±0.04 


0.57 ±0.03 


0.75 ±0.02 


0.61 ±0.03 


0.64 ± 0.03 


0.60 ±0.02 


Ha 


4.2 ±0.1 


3.9 ±0.1 


3.10 ±0.07 


2.9 ±0.1 


3.1 ± 0.1 


3.10 ±0.08 


[N II] A6583 


6.9 ±0.2 


5.7 ±0.2 


3.49 ± 0.07 


2.8 ±0.1 


3.4 ±0.1 


2.70 ±0.07 


[S II] A6717 


2.61 ±0.08 


1.74 ±0.05 


1.61 ±0.04 


0.85 ±0.04 


1.47 ±0.05 


0.89 ± 0.03 


[S II] A6731 


1.93 ±0.06 


1.19 ±0.04 


0.98 ±0.02 


0.58 ±0.03 


0.88 ± 0.03 


0.58 ±0.02 



a Coordinates of the aperture centers are given as offsets from the galaxy nucleus in arcseconds. 
b H/3 line flux in units of 10 -15 ergs -1 cm -2 . 




4000 5000 6000 7000 8000 

A (A) 



Fig. 1. — Integrated central spectra of A 3581 (top), 2A 0335+096 (center), and Sersic 159-03 
(bottom), summed over 5"-diameter apertures centered at the peak of the BCG continuum emission 
in each data cube. The spectra are corrected for the level of extinction indicated in each panel, 
derived as described in Section 13.21 
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Fig. 2. — Integrated and continuum-subtracted spectra from the cores of the BCGs; from top 
to bottom: A 3581, 2A 0335+096, and Sersic 159-03, showing spectral regions around the H/3 and 
[O III] lines in the blue data cube (left) and Ha and [N II] lines in the red cube (right). In each panel 
the upper solid line is the integrated galaxy spectrum, the dotted line is the continuum template 
spectrum, scaled to match the continuum within the displayed region of the spectrum, and the lower 
line is the continuum-subtracted result. Each spectrum is extracted from a 4"-diameter aperture 
at the galaxy nucleus. 



- 30 - 




AX (arc sec.) AX (arc sec.) AX (arc sec.) 



Fig. 3. — Left: A 3581 continuum image, centered at 8000 A in the galaxy rest frame. A foreground 
star is seen in the southwest corner of the image. Center: Map of the [N II] A6583 line emission 
surrounding the BCG. The cross mar ker indicates the locati on of the continuum peak. Right: 
Contours of 1.4 GHz radio emission from I Johnstone et al.l (|1998l ). overlaid on the greyscale emission- 
line image from our data. The radio beam size is shown in the bottom-left corner of the panel. 
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Fig. 4. — Velocity channel maps across the [N II]A6583 emission line. The velocity range in kms -1 
is shown at the lower left of each panel. The cross mark er indicates the positio n of the galaxy 
nucleus. Contours of the 1.4 GHz emission radio map from I Johnstone et al.l (|1998l ) are overlaid on 
several panels (black contours). 
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Fig. 5. — Left: 2A 0335+096 continuum image, centered at 8000 A in the galaxy rest frame. Center: 
Map of [N Hi line emi ssion surrounding the BCG. Right: Contours of 1.4 GHz radio emission from 



Donahue et al.l (|2007l ) are overlaid on the greyscale emission-line image from our data. In each 



panel, the BCG continuum peak is indicated by the black '+' marker, and the position of the 
companion galaxy by the ' x ' marker. 
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Fig. 6. — Properties of single Gaussian profile fits to the [N II] A6583 emission line for 2A 0335+096. 
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Fig. 7. — Velocity channel maps across the [N II] A6583 emission line detected from 2A 0335+096. 
The velocity range in k ms" 1 is indicated at th e lower left of each panel. Contours of the 1.4 GHz 
emission radio map from lDonahue et al.l (|2007l ) are overlaid on several panels (black contours). 
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Fig. 8. — Left: Continuum image of Sersic 159-03 centered at 8000 A in the rest frame of the galaxy. 
Cent er: Map of the [N 11 ] emission surrounding the BCG. Right: Contours of 5 GHz radio emission 
from iBirzan et al.l (|2008l ) overlaid on the emission-line image from our data. The continuum peak 
is indicated by the black cross marker. 
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Fig. 9. — Contours of the radio maps measured by iBirzan et al.1 (120041 ) from observations of Ser- 
sic 159-03 at 327 MHz (top left), 1.4 GHz (top right), 5 GHz (lower left), and 8 GHz (lower right), 
overlaid on a smoothed image of the [N II] emission from our data. The relevant radio beam sizes 
are indicated by the grey ellipses to the lower right of each map. 
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Fig. 10. — Maps of the line properties derived from single Gaussian profile fits to the [N II] A6583 
emission line, and, at lower-left, the [N II] /Ha line flux ratio for Sersic 159-03. 
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Fig. 11. — Velocity channel maps across the [N II] A6583 emission line in Sersic 159-03. The velocity 
range in kms -1 is indicated at the lower left of each panel. 



